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(54) Linear space-time muttluser detector 

(57) Linear space-time nujltiuser detectors are 
described that combine an'ay processing, rake detec- 
tion, and multiuser detection and that satisfy a minimum 
mean-squared error criteria. Both embodiments can be 
inplemented adaptively so as to account for unknown 
sources of interference. The adaptive implementation of 
the first embodiment of a linear space-time detector 
requires explicit estimates of the channel coefficients. 
On the other hand, the adaptive implementation of the 
second embodiment does not require these estimates 
since they can be adaptively obtained. Both of these 
detectors provide significant performance advantages 
over single-user space-time detectors and over multi- 
user detectors whch do not account for CDMA interfer- 
ence from unknown sources. 
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Description 

Reld of the Invention 

5 [0001] This invention relates generally to direct sequence (DS) code division multiplex access (CDMA) systems sub- 
jected to muKipath fading and multiple access interference (MAI), and more particularly to signal detection techniques 
for use in such systems. 

Background o1 the Invention 

10 

[0002] In wireless CDMA systems such as those proposed for next generation mobile wireless standards, the major 
impediments at the physical layer are muttipath fading and multiple access interference caused by co-channel users 
which are not orthogonal to the desired user. The rake receiver, described in R. Price and P. E. Green. "A Communica- 
tion Technique for Multipath Channels/ Proceedings of the IRE. Vol. 46, pp. 555-570, March 1958, attempts to combat 
75 multipath fading by coherently combining resolvable multipath replicas of the desired signal. Multiuser detection 
descrik>ed in S. Verdu, "Multiuser DetecthnT Cambridge University Press. New York. 1998, addresses the problem of 
MAI by actively accounting for its presence when detecting the desired user. 

[0003] More recently, there has been growing interest in using array processing for further improving receiver perform- 
ance. These techniques have focused on using multiple arrtennas at the base station receiver to provide antenna gain 

20 arxi/or diversity gain and allow the possibility of spatial processing. By combining these space-domain techniques with 
time domain techniques like rake detection and multiuser detection, the resulting space-time detectors show promise 
of improving the capacity of CDMA systems as compared to traditional time-domain-only detectors. See, for example. 
A. Paulraj arxi C Papadias. "Space-Time Processing for Wireless Communications." IEEE Signal Processing Maga- 
zine. Vol. 14, No. 6. pp. 49-83, Nov. 1997. The first generation of space-time CDMA detectors used array processing 

2S with either rake detection or multi-user detection. See. respectively, A. Naguib and A. Paulraj, "Performance of Wireless 
CDMA with Af-ary Orthogonal Modulation arxl Cell Site Antenna Arrays." IEEE Journal on Selected Areas in Commu- 
nications, Vol. 14, No. 9. pp. 1770-1783, Dec. 1995 orS. Miller and S. Schwaru. "Integrated Spatial-Temporal Detectors 
for Asynchronous Gaussian Multiple-Access Channels." IEEE Transactions on Communications, Vol. 43. No. 2/3/4, pp. 
396-411. Feb./Mar7Apr. 1995. Later space-time CDMA detectors combined all three processing techniques. See: H, 

30 Huang; S. Schwartz. S. Verdu. "Combined Multipath and Spatial Resolution for Multiuser Detection: Potentials and 
Problems." Proceedings of the IEEE International Symposium on Information Theory, p. 380, 1995; or M. Nagatauka 
and R Kohno. "A Spatially and Temporally Optimal Multi-User Receiver Using an Array Antenna for DS/CDMA" lEICE 
Transactions on Communications, Vol. E78-B. No. 11. pp. 1489-1497. Nov. 1995. 

[0004] While the foregoing systems operate satisfactorily, improvements can be made. In particular to the space-time 
35 detectors which combine all three processing techniques. While the detector in the Nagatsuka and Kohno paper is opti- 
mum in the maximum likelihood sense, its computational complexity is exponential with respect to the number of users. 
Hence it is too complex to implement for practical systems. In the paper by Huang, Schwartz and VerdO. a tradeoff 
between performance and complexity is made, but tNs detector was not implemented adaptively since it used a zero- 
forcing criteria. Adaptive Implenerrtations allow receivers to account for unknown sources of interference thus Improv- 
40 ing the detector performance and increasing the system capacity. For exanple. a base station receiver could account 
for interference from adjacent cells or from an embedded microceQ, while a handset receiver could account for interfer- 
ence from signals it is not expfidtly demodulating. 

Summary of the Invention 

[000$] This Invention detects DS-CDMA signals utili2lng a rake receiver, array processing, arid multiuser detection. 
When combined with array processing, the rake receiver is often caDed a space-time rake receiver. The invention uses 
a minimum mean-squared error (MMSE) aiterion in the multiuser detector. This criterion allows for relatively simple 
implementation in the form of a linear detector arxl also allows tor adaptive implementation. Adaptive implementations 

so are useful in practical situatiors where there is limited knowledge of the various received signals. As (tiscussed above, 
both uplink and downlink capacity can be improved using adaptive detectors which account for unknown interference. 
[0006] For pulse amplitude modulated (PAM) data signals, two emtxxJiments of the invention provide options for trad- 
ing off between performance and adaptive implementation complexity. The first embedment p&forms better if there is 
perfect knowledge of the users' signal parameters. On the other hand, it requires wore explicit channel information for 

55 adaptive implementation. Suppose the receiver consists of 1 antennas arxl the received signal consists of /C DS- 
CDMA data signals, each with L detayedMeighted multipath replicas. At each antenna, a bank of filters is matched to 
the KL spreading codes with their multipath timing delays. The fitter outputs are weighted according to the complex con- 
ji^te of the estinriated channel (muttipath and array) parameters and combined to form a K-vector wfiere each oom- 
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ponent corresponds to one of the K codes. The real part of each component is taken. Estimates of the nrujHipath delays 
and channel parameters may be obtained, for example, from a training or pilot signal. The conventional space-time rake 
receiver would pass each vector component to a decision device for estimating the corresponding user's PAM data 
symbol. However, unless the users* signals are orthogonal in the space-time domain, the vector components are con- 

5 tarrenated with multiaccess interference from the other users. This invention proposes the use of a linear combiner prior 
to the decision devices to suppress the multiaccess interference. This linear combiner is represented by a real K-by-K 
matrix W/^ which multiplies the real K-vector. The matrix minimizes the mean squared error between its product with 
the real K-vector and the K-vector of data symt>ols. Because the matrix uses the minimum mean-squared 
error(MMSE)criterion, well-known adaptive algorithms can be used to adaptively obtain it. Each component of the final 

70 K-vector output corresponds to one of the K codes and is passed to either a decision slicer or a decoder for further 
processing. 

[0007] The second eir^odiment has marginally inferior performance when the channel parameters are known exeictly. 
However, under practical conditions of channel mismatch, the second embodiment often performs better. In terms of 
adaptive implementation, this errt^odiment requires less information (rt does not require explicit channel estimates) tnJt 

IS may be slower to adapt. As vM\ the first embodiment, the front end of the second embodiment consists again of a bank 
of KL matched fBters for each of the P antennas followed by weighting and combining. However, unlike the first embod- 
iment, the order of the real operator and linear combiner(K-by-K matrix multiplication) are exchanged. In doing so. the 
weighting, combining, and linear combiner can be represented by a single complex K-by-KLP matrix Wq, The real parts 
of the components in the resulting K-vector are passed to either a decision device or a decoder for further processing. 

20 In the adaptive implementation, adaptive algorithms can be used to obtain the matrix W^: hence channel estimates are 
not explicitly required. 

[0008] For quadrature PAfS/l (QAM) data signals, real operators are not required in the detector since the signal con- 
stellaton is two-dimensional. In this case, a variation on the secorK) emtxxiiment would be used which does not use the 
real operators but follows the liner combiner directly with the appropriate decision device or decoder. 

25 [0009] While the I^MSE detection techniques are powerful in and of themselves, their performance can be further 
enhanced by using them with other multiuser detection techniques such as interfererKe cancellation where interfer- 
ence, instead of being projected away as is done with linear multiuser detectors, is explicitly subtracted from the 
received signal. Interference cancellation can occur prior to or after the MMSE linear combiner. 
[0010] Both embodiments of the Mfy/ISE detector can be generalized for practical systems which utilize, in addition to 

30 the data bearing channels, supplemental channels which may serve as pilot signals. Also, the IVIMSE detectors can 
operate in systems where signals are transmitted with different spreading factors. 

[0011] A CDMA signal detector in accordance with the invention can provide significantly improved performance rel- 
ative to a conventional signal detector which does not use multiuser detection or which does use multiuser detection 
but does not adaptively mitigate unknown interference. 

35 

Brief Description of the Drawings 

[0012] The foregoing and other features of the invention may become rnore apparent when the ensuing description 
is read together with the drawings, in which: 

40 

Figure 1 shows a prior art single-user space-time rake receiver employing a multi-element antenna array, a bank of 
N'Ctup filters at each antenna, and a channel weighter and combiner; 
Figure 2 shows the generic embodiment of the inventton; 

Figure 3 shows a first embodiment (Detector A) of the multiuser space-time MfS^SE detector, 
45 Figure 4 shows a second embodiment (Detector B) of the nrubuser space-time MMSE detector. 
Rgure 5 shows an adaptive implementation of the second embocGment; 

Figure 6 shows interference canceOation located prior to the linear combiner of the first embodiment; 
Figure 7 shows interference canceOation located prior to the linear combiner of the second embocfiment; and 
Figure 8 shows interference cancellation foDowing the linear comt»n&^ of the multiuser space-time MMSE detector 
so in accordance with the invention. 

General Description 

[0013] Conskier a K user system where the Mh user (/c » 1 ... K) modulates its data sequence bi^f) with an N-chfp 
55 spreading sequence Sj^Q. For Detector A the data sequerKe is pulse amplitude nrvxjulated (PAM); for Detector B. the 
data sequence can be other PAM or quadrature PAM (QAM).. The transmitted signal uridergoes frequency-selective 
fading in the channel and arrives at the receiver as L lime-resolvable multipath components with complex fading chan- 
nel coeffidents q^/(0 "C*;i.W* Assume that the receiver is a P-element linear array. If each resdvabie mult^th com- 
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ponent arrives as a planar wavefrorrt with angle G f^f{f) with respect to a linear array and if the array spacing is 
sufficiently close (e.g., }J2) so that there is perfect correlation annong array elements for a given wavefront then the 
phase offset of the pth element, with respect to the first is /7^^(0=e3cp(ii(p-1)sine ^;(0) - The received signal at the 
pth antenna for given symbol period (and ignoring intersymbol interfa-ence) is: 



lr«1 /»1 



(1) 



10 



where Af^ is the amplitude for user k. t/f / is the delay for the kth user's fth multipath. and np{t) is the additive Gaussian 
noise process which accounts for out-of-cell interference and background rroise. We now make the following assump- 
tions to simplify the analysis. 



(a) The signals are received bit synchronously. 
75 (b) The time spread is small compared to the symbol period so that intersymbol interference can be ignored.: and 
(c) The phase offsets and channel coefficients are constant over a symbol period. 

The first assumption will be dropped later. Under these assumptions, the chip matched-filter output for the received sig- 
nal in equation (1) is a complex A/-vector 

20 . .. 

25 

where s^j is the chip matched-filter N-vector corresponding to Sff(t -r^f /) and Hp. is the conplex N-vector corresponding 
to the Gaussian noise. The spreading codes are normalized to have unit energy: 



30 



The spreading codes are assumed to be random. However, adaptive implementations of these fVir^SE detectors require 
35 the use of short spreading codes which repeat after a few symbol periods. The following notation will be used hereinaf- 
ter: 



40 h'^s c/>ap(hi^V"hi,i.. 

C= tf/ap(Ci...CK) 
R 



complex P-vector of an-ay coefficients 

complex PxKL array matrix 

complex KLP x KL array matrix 

complex L -vector of channel coefficts 

complex KLxK channel matrix 

complex KL x KL correlation matrix defined by: 



45 



so A = diag{fii^..Af() 



ss 



real K X /C amplitude matrix 
real K-vector of data bits 
u X u identity matrix 
u-vector of ones 
uxu matrix of zeros 
u-vector of zeros 



The corresponding estimated values for the array and channel coefficients will be r^resented with ttie symbol over 
the symbol for the value. The nodse vector is a zero-mean complex Gaussian vector whose di^ritxjtion can be written 
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in terms of its (component-wise) real and imaginary components: 



Rc(n) 
Ini(n) 




where we define the real and imaginary operations for matrices and vectors as 



10 



Rfi(X)=(X+3^)/2 and lii<jg=(X-X')/2 



IS where * denotes the complex conjugate. Hence Re{n) and lm(n) are zero-mean Gaussian random vectors whose com- 
ponents have variance and are mutually uncorrelated. 

Conventional (Space-Time Rake) Receiver, (Rgure 1) 

20 [0014] In the context of this array-multipath channel, the prior art single user detector is a correlator matched to the 
desired user's composite array-multipath-spreading code signal. This detector does not account for the presence of 
interterers; however, it is the maximum-likelihood detector if there are no interferers or if they are orthogonal to the 
desired user in array-code space. As shown in Figure 1. the detector consists of a bank of correlators 10 at each 
antenna, matched to the KL multipath spreading codes s-j v- Si ^/....SK.v .SK.i- The notation < /; )in the left-nrrast boxes 

2£ Of Rgure 1 indicates taking the inner product between s^,] (the icth user's spreading code corresponding to the rth mul- 
tipath delay) and the input vector to the box. The dot in the notation (S/(^/;> represents the input. The timing estimates 
are obtained i^ing a separate timing estimation algorithm. The output of the upper left-hand box (which corresponds to 
the first component of the vector 2^) is 



where the H superscript signifies the Hermitian transpose (take the complex conjugate of each element, then take the 
36 transpose of the resulting vector) of a complex vector. The matched fSter outputs for at least a subset of the multiple 
signals, a plurality of the multipath components of these signals, and a plurality of receiver antennas are indicated in the 
drawing at (z^ ...Zp) 

[001 5] At reference number 1 1 . the correlator outputs are weighted by the complex conjugate of their corresporKfing 
channel estimate. I^ore specifically, the conetator output at the pth antenna for the /th multipath of user k is weighted 

40 by the complex conjugate of the estimate of the corresponding channel (array/multipath) coefficient h kjpC^j. These 
estimates are obtained using a separate channel estimation algorithm. At each antenna, the L components for the /cth 
user are added (l 2). and then the resulting P components for the Mh user are added 1 (>c « 1 ... K). Each component 
is then passed to either a decision device or a decoder. For PAM data signals, the decision device tor the Mh user (14^^. 
ks 1 ... K) gives as its output the symbol which is closest in the Euclidean distance smse to the real part of a /cth vector 

45 component. For QAM data signals, the real operator is not required. 
[0016] From Rgure 1, the K/J'-vector 



30 




SO 




at the output of the matched fDters can be written as: 



£5 



(2) 
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where we define the KLP x KLP matrix 

R=R0I^ 

and the KLP x KL matrix . ^ denotes the Kronecker product operation between two matrices, and the complex noise 
vector is a complex Gaussian random vector with distribution 



Un')J [-In^) Ri^)JJ 



100171 Using the fact that 

(H''c)''RH'^C = e'[(H''H)oR]c 

(where . denotes the component-wise product between two same-sized matrices), the K vector at the output of the 
channel conrd3iners and the input to the decision devices 14^-14|^ can be concisely written as: 

Mya».}=Re{(H''cfz}=MAb+D_ (3) 

where H denotes the complex conjugate transpose, ° represents the component-wise product of two same-sized matri- 
ces, and 

MsRe{c"[(H"H)oR]c} 



The noise vector is purely real and has distrObution: 



where 



M=C*[(H''H).R]e. 

H the data is BPSK modulated, the tM decision for user k is eirrply the hard Emit of the kth component of 
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5 



The corresponding bit error rate P^"^ of the conventional rake recever is: 



10 




(4) 



where y^(k,k) »s the (KK)th element of matrix X. The performance comparisons among the various detectors are given in 
terms of BPSK modulated data, but they can also be done for any QAM mocftjiated data in general. 

20 Generic Space-Time Linear Multiuser Rake Receiver (Figure 2) 

[0018] The prior art space-time rake receiver described above is a single user receiver in the sense that for a given 
user, the demodulation uses information only from that user. Since it does not account for the presence of interference 
from other users, its performance suffers. Rgure 2 shows a generic multiuser space-time rake receiver. It accounts for 

25 the multiaccess interference using a linear combiner. Two versions of this detector are given. The first, known as detec- 
tor A. uses the real operators before the linear combiner. In this case, reference 20 in Figure 2 consists of a bank of real 
operators followed by a linear combiner, and each of the components 14^.-1 V is a slicer which determines the closest 
estimated symbol to its input. The second, detector B, uses a linear corrdsiner represented by a complex matrix followed 
by the real operators. In this case, reference 20 consists only of the linear combiner, and the decision devices consists 

30 Of a real operator and slicer. Despite the similarities, the fact that the linear combiners for detectors A and B results in 
different performances and adaptive implementations. We now describe the two detector embodiments more thor- 
oughly. 

Space-Time Linear IMultiuser Detector A (Rgure 3] 

35 

[0019] A set of suffiderrt statistics for the K users is given in equation (3) as the K- vector Re{yjponv}- The objective of 
the linear MMSE detector is to apply a fC x AC linear combiner to this vector such that the mean-squared error between 
its resultir^ vector and the data vector b is minimized. In other words find the real KxK matrix W/^ such that 



This detector in essence strikes the prefer balance in oorrbatlng both residual Gaussian noise, manifested tlvough 
45 n^"^ . and multple access interference, manifested through the off-diagonal terms of 



40 




(5) 




SO 



It can be shown that the solution Is: 



ss 
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w 



.£[b(R«{y-}f]{4Re{r"}R.''{y-}]}" 



= AM(mA W +a^M)"^ 



[0Q20] Figure 3 shows a Wock diagram of detector-A in accordance with the invention. It contains P antennas and a 
bank of matched fitters 10. The box 31 represents a linear combiner d^ined by the K-by-K matrix Wyj given by equation 
(6). The matrix can be put Into words as follows: it is the product of 3 terms: (a) the diagonal anriplttude matrix (A); (b) 
the equivalent array/channel/code correlation matrix (M); and the matrix inverse of the matrix 



vtfhich consists of the product of the array/channel/code correlation matrix (M). the square of the amplitude matrix (A^). 
and the Hermitian conjugate of the array/channel/code correlation matrix (M^ summed with the estimated array/chan- 
nel/code correlation matrix 



which is weighted by the background noise variance (o^. This matrix minimizes the mean-squared enor criterion given 
in equation (5). and because of the convexity of the mean-squared error expression in (5). the MMSE solution (6) rep- 
resents a global minimum. A bank of decision slicers (32i - 32k) follows the linear combiner. 

[0021] An attractive feature of the MMSE detector is that it can be implemented adaptively using well-known adaptive 
algorithms (33) such as least-mean-squares or recursive least-squares. Using a training signed (for example, the user 
data bits), these techniques can be used to adaptively obtain W^. Adaptive implementation Is an option If direct calcu- 
lation of Is deemed too corrplex. Note that the adaptive inplementation requires knowledge of the K users* spread- 
ing codes, multipath delays, arxj channel parameters. 
[0022] Assuming that the array and channel estimates are exact 



and another attractive feature is that MMSE detector approaches the decorrelating (zetD-forcing) detector {vp to a sca- 
lar facto-) as the residual Gaussian noise approaches zero: 



(MAW+(7^M) 




C=C,H = H. 



we can rewrite the expression for Wy^ in (6) as: 






I n other words, from (3) we have that 
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= b 



10 and the multiple access interference is forced to be zera In the implementation of the decon^elating detector A (DD-A). 
we replace Wy, given by equation (6) with 



IS 



20 



W^=A-*Re-'{R^. 



We will assume that C^[(H^H)i^lC is positive definite and hence invertible; and it follows that Re{C"((H"H)J?]C) is 
also invertible. 

[0023] Let us now calculate the bit error rate of this first f^MSE detector (MMSE A). The K-vector at the slicer inputs is 



25 where 



30 



Hence the bit error rate for the klh user using the MMSE-A detector is: 



35 



(8) 



40 



46 



Letting 



W, = A-'(Re{c*[(H*H)cR]c})"' 



and assuming perfect channel estimates, the bit error rate for DD-A, is: 



55 
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[0024] One drawback of this MMSE detector is that the estimation of the channel and array coefficients cannot be 
incorporated into the adaptive algorithm for obtaining W^. The estimates must be explicitly obtained using some other 
means (like a training or pilot signal). This leads us to ask. "is it possible to design an adaptive MMSE which does not 
require prior knowledge of the array and channel parameters?" The answer is a"yes." as we will see in the following sec- 
tion. 

Space-Time Linear Multiuser Detector B (Rgure 4) 

[0025] The goal of this subsection is to derive an MMSE detector whic^ can be implemented adaptively and which 
does not require explicit array and channel estimates. (It will, however, still require timing estimates for the muttipath 
delays.) Figure 4 shows the basic structure of Detector B which exchanges the places o1 the real operators and linear 
combiner of Detector A. Following this exchange, the resulting cascade of the array combining, multipath combining, 
and linear combiner can be collected into a single complex K x KLP matrix multiplication called W^. as shown in Figure 
5. Specifically, we define this matrix using the MMSE criterion as: 

W, = argimii{£[||Wz-b|' 

and it follows from (4) that 

W, = £(!«-)[£(«-)]-' 

^AG^RlRGA'C^R+iT^R]"*. (10) 
= AG''[RGAH5''+<r'I^f 

where 

is the ray/channd matrix. Figure 5 shows a block diagram of the MMSE-B detector in accordance with the invention, ft 
contains P antennas and a t>ank of matched filters 10. The box 50 r^resents a linear combiner defined by the K-by- 
KLP matrix given by equation (10). This matrix mininnizes the ntean-squared error criterion given above by 

W, = argmm|£[||W2-bf|. 

The last line of equation (10) states that the matrix is the product of 3 terms: (a) the diagonal amplitude matrix (A); (b) 
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the array/channel combiner matrix H^C which is equivalent to the array/channel combiner in the conventional detector 
but using the actual parameters and not their estimates; and (c) the matrix inverse of the matrix 



The matrix inverse consists of the product of the correlation matrix 

(R) 



defined following equation (2). the ray/channel combiner matrix G, the square of the anplitude matrix (A^), and the Her- 
mitian conjugate of the array/channel combiner matrix summed with the KLP-by-KLP identity matrix weighted by the 
background noise variance (cf^. ff the data is PAM.. a bank of real operators (30-j - 30k) decision slicers {32-| - 32k) 
follows the linear combiner. Otherwise, if the data is QAM. only the complex decision slicers (14^ -14^) are required. For 
PAM data. Detectors A and B, while similar in structure, are not equivalent due to the different placement of the real 
operator. 

[0026] To calculate the bit en^or rate of Detector B for BPSK data, the K-vectbr at the slicer inputs is 



where 



Hence the b*t error rate for the kih user using the MMSE-B detector is: 



^c(tVjRH''c)u>} 



[0027] By concatenating the P received /V-vectors into a single A/P-vector 

r=(,f^.,.^]'. 

we can combine the operations of the matched fitter bank with the linear combiner to create a single K-ty-NP 
matrix which multiplies r. As shown in Figure 6, the detector structure for the kth user can be reduced to the /cth row of 
this K-by-NP matrix 60 followed the real operator 62 and decision sGcer 63. Hence the MMSE-B detector can be 
implemented with a simple tap-weight filter ardirtecture. Its corresponding decorrelating detector (descra^ed below), the 
MMSE*A detector arxJ the DD-A detector can likewise be implem^rted with tiiis architecture. Like the MMSE-A detec- 
tor, this version of the MMSE-B detector can be implemented adaptively using well-known adaptive edgorithms (61). 
However. MMSE-B has the advantage of not requiring the array and channel estimates. The tradeoff ie that the adap- 
tation may be slower t>ecau&e there are more taps to acfjust 

[0Q28] While direct analytic comparisons between the two detectors is difficult, we will show that under perfect chan- 
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nel and array estimates, the performance of MMSE-A is uniformly superior to that of MMSE-B as the noise floor 
approaches zero. First, we show that the corresponding decorrelating detector for MMSE-B, DD-B. simply replaces 
with 

(Recall that we have already made the assumption that 

10 

(h^c)''rh''c = c"[(h''h)or]c 

75 

is invertible.) 
Proposition 1 : 
20 [0029] 

liJA(H^C)''[RH''CA'(H^C)" +0^1^]"'!^ 



Proof: Multiplying tx}th sides by 

30 

35 we have that 

A(H''cY = limA- j(H^C)*ra^c]"'(H*C)''[RH''CA'(H"C)' + a^l^^ 

= a(H^C)* + ljm(cT*A-'[(H''C)*to''c|"*(H''C)'') 
»A(H''cf. 

Hence, from the fact that 

{H^C)''fiH'*C=C''[(H''H)oR]c 

and the fact that the array and muttipath combiners in MMSE-A can be represented by (H^D)^. the DD-B is the bank 
of matched filters, followed by the combiners, followed by the decorrelator (C^[(H'^H}i)]C)'^ . followed by the real oper- 
55 atbrs. and followed by the slioers. The DD-A is nearly identical exc^ that the order of the deconelatdr and the real 
operators are switched, and the decorrelator is (Re{C^[(H^H)* R]C})'^ 
By replacing with 



RH^CA'(H^cf 
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5 



in (1 1), the bit error rate of DD-B is: 



10 



A, 



(12) 



IS 



aJlu[{c"[{H''H).R]c}-] 



^0 [0030] We would now like to show that DD-A has a uniformly lower bit error rate than DD-B. 
Proposition 2: P^^^ (a) <. P^^^ (a) 

[0031] Proof (for which the inventors wish to thank Emre Telatar of Lucent Technologies* Math Sciences Center for 
25 his timely and elegant insight): 
[0032] Define 



which is Hernnitian and assumed to be positive definite. It follows that [Re(M)]'^ is also positive definite, and it is suffi- 
cient to show that the kth diagonal element of [ReCIDf)]'^ is greater than or equal to the kth diagonal element of IR'''. Vt 
V (Note that since 19 is Hermitian, its inverse is also Hermitian, hence the diagonal elements of the inverse are real.) 
35 Chose y to be real and set 



30 



MsC"[(H''H)oR]c, 




40 



and note that X is real. Thus, 



4S 




ami also 



so 




Using a generalization of Bergstrom's inequaGty whi^ states that for any real x, real y, and positive def inite^n. 



55 
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we have that 



10 



25 



30 



y'[Re(Ri)]-yy'(JS)-'ys(y'[RKIB)rV)' 
y'(ia)"'y^y'[Re(ia)]''y. 



and the result follcfws by using the kth unit vector for y. Desp'rte this relationship, it turns out that the perforniance of DD- 
A and DD-B are actually quite similar. More importantly, it has been shown that the performance of MMSE-B can be 
better than MMSE-A for the practical case of inexact array and channel estimate. 

75 Extensions of MMSE Detectors to the Asynchronous Multirate Case 

[0033] So far, we have assumed bit synchronous operation of these linear multiuser detectors. We now show how 
their operation can be extended to the bit asynchronous case. First, we consider the operation off a decorrelating detec- 
tor in a synchronous singte-rate environment. Similar reasoning applies to MMSE detectors; however, for pedagogical 
20 purposes, we focus on the decorrelating detector. Suppose that there are two users with respective spreading codes 
s^(0 and 52(0* Table A below shows the relationship between the two codes in time, with the thin vertical lines repre- 
senting symtxjl boundaries. The received signal Is r[t)^A^s^{t)b^ + ^252(0^2 '^(0 where is the Mh user's 
anplitude, bf^ is the /cth user*s data bit. and n{f) is the additive Gaussian noise. 



• • 


*.(/) 


9 « 


• • 




• « 



35 



Table * A: 2 user synchronous single-rate system 



40 



46 
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SS 
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Table - B: 2 user asynchronous single-rate system 

10 
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• • • 


• • • 








-27^ 


• • • 



Table - C: 2 user asynchronous multi-rate system 



20 



[0034] H user 1 is the desired user, its decorrelating detector lor a given symbol interval is the nnatched litter for code 
Si(0 projected into the null-space of S2(0- For the asynchronous single-rate case, shown In Table - B, the decorrelating 
2£ detector for user 1 is projected into the null-space spanned by the linear combinations of user 2 s codes which overlap 
during that symbol interval. H we let s « (0 be the portion of SgCO which overlaps with (0 and let (f- 7) be the portion 
of S2(^-'0 which overlaps with Si(f). then, assuming BPSK data modulation, the null space is spanned by 

s^(0+4<'-"n ands^{f)-S2(^r). 

30 ... 

This idea can be extended to the asynchronous multi-rate case shown in Table - C Using similar function definitior^. 
the decorrelator for user 1 lies in the null space spanned by 



35 



40 



s^(t)^s^{t-Ty,s^{t'2n s^(0+S2(f-r)-s^{/-2T). 

As the disparity between data rates increases, the eHectiveness of a linear multiuser detector decreases since it will be 
constrained to a more restrictive subspace. The reasoning above can be applied to MMSE detectors by simMarly com- 
bining the appropriate subspaces of the interfering waveforms. 

MMSE Detectors and Interference Cancellation 

45 [00351 The MMSE techniques described here can be used inconjunction with a non-linear multiuser detection tech- 
nique known as interferwce cancellation. \Whereas the linear multiuser detectors rely on subspace projections to miti- 
gate interference, the interference cancellers directly subtract off the interference. Two types of interference cancellation 
can be used: pre-con*iner, vtfhere the cancellation occirs before the linear combiner, and post-combiner, where it 
occurs after. Pre-con*iner cancellation is used for canceling interference which is not accounted for by the linear corn- 
so biner. For exarrple. we may choose to cancel a high-powered, high-data rate signal from within the cell instead of 
accounting for it with the linear combiner. Post-combiner cancellation on tiie other hand is used for refining the 6yrT*>ol 
estimates made by the linear conribiner. Using preliminary symbol estimates from the linear multiuser detector, the sig- 
nals are reconstructed and subtracted from the received signals to form enhanced signals which are used for a second 
stage of symbol estimatioa In both interference cancellation techruques. capacity can be potentially increased by sup- 

55 pressing the irrterference. ^ * u 

[00361 Figure 7 shows a block diagram of a pre-combiner interference canceller in accordance with Detector A shown 
in Rgure 3. At the Mh input to the finear combiner W^. there is interference due to users / « 1 ...K. j ^ k. which wiQ be 
accounted for with W/^. However, there may also be othenwise unaccounted interference, for example, high-powered 
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interference from within the cell. This interference could be subtracted from the ir^ut to Wy^ if the spreading codes, mul- 
tipath delays, data bits, and channel parameters at all antennas are known. In this case, the baseband signal for the 
interference could be recor^tructed and its contribution to the k\h input of W/^ could be estinnated (by passing the recon- 
structed signal through the processing c^ain for user Ac) and subtracted at the poirrt designated 70j^ (Ac = 1 ,.. K) Contri- 

£ butions from other interference for this or other inputs can be similarly computed and siArtracted. 

[0037] Figure 8 shows a block diagram of a pre-combiner interference canceller in accordance with Detector B from 
Figure 4. Using the same method as desaibed above for the Detector A pre-combiner interference canceller, the inter- 
ference which is not accounted for by can be subtracted from each of its inputs at the point designated 80. 
[0038] Figure 9 shows a block diagram of a post-combiner interference cancellation receiver in accordance with the 

10 invention. Ths receiver includes P antennas and has three stages, designated 90. 91 and 92. The first stage 90 
receives baseband signals fp, and makes preliminary symbol estimates for all K users with either a MMSE-A detector 
of Figure 1 or a MI^SE-B detector of Figure 2. The second stage 91 uses the preliminary symbol estimates from the 
first stage 90 and knowledge of the K users' spreading codes, delays, and channel parameters to reconstruct the base- 
band received signals for each of the users. The second stage 91 subtracts the multiaccess interference with respect 

15 to a desired user k from the received signals r^...rp to form the following enhanced received signals for user Ac. 



[0039] Note that the multipath interference with respect to user k can also be removed. We will not describe this option 
in detail because its relative gains will generally be negligible unless the number of users K is very small. Ideally, if the 
2& preliminary symbol estimates and the channel estimates are perfect, there will be no multiaccess interfererKe in the 
enhanced received signals. However, this will generally not be the case in practice. In the third stage 82. the enhance 
received signals are processed with a conventional space- time receiver shown in Figure 1 to generate the final symt>ol 
estimates b J,'^) for the K users. 

[0040] The post-combiner interference cancellation procedure described above can be continued for multiple itera* 
3D tions by repeating stages 91 and 92 indefinitely. Performance gains due to repeated iterations diminish after the first two 
or three iterations. To reduce the complexity, interference cancellation of both types can be performed on a subset of 
users. 

[0041] What has been described is deemed to be illustrative of the principles of the invention. In particular, it should 
be noted that the apparatus and methods of the invention may be implemented through various technologies, for exam- 
35 pie, through the use of large scale integrated circuitry, application specific integrated circuits and/or stored program 
general purpose or special purpose computers or microprocessors using any of a variety of computer- readable media. 

Claims 

40 1. A method for detecting the data signals modulated respectively by multiple direct-sequence spread spectrum sig- 
nals, comprising the steps of 

a) generating nrtatched filter outputs for at least a subset of the nrudtiple signals, a plurality of the muHipath com- 
ponents of these signals, and a plurality of receiver antennas; 
45 b) de\/elpping coherent channel estinnates: 

c) weighting and combining the matched fitter outputs using the coherent channel estimates; and 

d) estimating the respective data symbols for at least a subset of the spread spectrum signals while suppress- 
ing multiple access interference. 

so 2. The method of claim 1 wherein said coherent channel estimates are developed from said matched filter outputs. 

3. The method of claim. 1 wherein sakJ coherent channel estimates are developed from an auxi&ary signal. 

4. The method of claim 1 wherein multiple access interfererice is suppressed with a iinear combiner and wherein a 
55 decision device optimally estimates the respective data symbols given the linear combiner outputs. 

5. The method of dalm 4 wherein the steps of generating matched filter output; of combining the matched filter out- 
puts, and of utilizing a linear comb'n^ are implemented using a modified tap-wdght f 3ter. 
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6. The method of claim 4 wherein the linear corr^iner employs a minimum mean-squared error criterion. 

7. The method of daim 6 wherein the nruntmum mean-squared error criterion designates a linear combiner for the Mh 
user which nrunimizes the expected squared error between the combiner output and the Mh user data syrr4)ol: 

5 ' ■ ■ 

8. The method of claim 7 wherein the linear combiner forms a matrix V/q from the complex block-Toeplitz correlation 
matrix 



10 




the diagonal amplitude matrix A, the complex ray/channel matrix G. and the background noise variance o^. 

9. The method of claim 7 wherein the linear corr^iner forms a matrix from the complex block-Toeplitz correlation 
i£ matrix 

ft. 

20 the diagonal amplitude matrix A, an estimate of the complex array/channel matrix 

Q 

25 and the background noise variance o^. 

10. The method of claim 4 wherein the linear combiner employs a zero*forcing error criterion. 

11. The method of claim 10 wherein the linear combiner forms a matrix from the complex biock-Toeplitz correlation 
30 matrix 

R, 

36 

the complex ray/channel matrix G, and the diagonal amplitude matrix A. 

12. The method of daim 10 wherein the linear combiner forrr^ a matrix Wg from the complex block-Toeplitz correlation 
matrix 

40 

ft. 

an estimate of the complex array/channel matrix 

46 ' 

G. 

and the diagonal amplitude matrix A. 

so 

13. The method of daim 4 wherein the data signals are pulse anplitude modulated, and the real components of each 
of the combined matched filter outputs are extracted prior to the linear combiner. 

14. The method of daim 13 wheran the linear combin&^ enrploys a minimum mean-squared error criterion. 

55 

15. The method off daim 14 wherein the minimum mearvsquared error criterion designates a Gnear connbiner for the 
klh user which minimizes the expected squared error between the comt>iner output and the /cth us& data symbol. 
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16. The method of claim 15 wherein the linear combiner forms a matrix W^, from the real equivalent array/channet/code 
correlation matrix Re(M], the diagonal amplitude matrix A. and the background noise variance o^. 

17. The method of daim 15 wherein the linear combiner iorws a matrbc Vi^ from the real part of an estimate of the 
£ equivalent array/channel/code conelation matrix 

Re(lVf). 

10 the diagonal amplitude matrix A, and the background noise variance o^. 

18. The method of claim 13 wherein the linear combiner employs a zero-forcing error criteria. 

19. The method of claim 18 wherein the linear combiner zero-forcing criterion includes forming a matrix Wy, from matrix 
T£ inverse of the real equivalent array/channel/code correlation matrix Re(M). 

20. The method of claim 18 wherein the linear combiner zero-forcing criterion includes forming a matrix W>; from the 
matrix inverse of the real part of an estimate of the equivalent anay/channe!/code conelation matrix 

20 / ^ \ 

Re(M). 

21. The method of daim 1 wherein the direct sequence spread-spectrum signals are spread with different spreeuJing 
gains. 

22. The method of claim 4 wherein, prior to employing the linear combiner, interference cancellation is performed on at 
30 least a subset of the combined matched fitter outputs. 

23. The method of claim 4 wherein the linear combiner is followed by interference cancellation on at least a subset of 
the linear combiner outputs. 

3£ 24. Apparatus for use in delecting at a receiver in a wireless system with multiple ti^ansmitted direct-sequence spread 
spedrum signals, comprising: 

a) means for generating matched filter outputs for at least a subset of the multiple signals, a plurality of the mul- 
tipath components of these signals, and a plurality of receiver antennas; 
^0 to) means for developing coherent channel estimates from the matched filter outputs: 

c) means for weighting and combining the matched fitter outputs using the coherent channel estimates; and 

d) means for estimating the respective data signals for at least a subset of the spread spectrum signals while 
suppressing multiple access interference. 

4S 25. A receiver comprising: 

a) space-time rake receiver having: 

i} a k>ank of correlators for resolving K code-spread data signals Si^....,5/c bearing respective data 

so b^,bz each with at rnosi L resolvable muttipath components with complex amplitude Cj^ (Ac « 

1,2....,K; / B 1.2....,/.), arriving at a P element antenna array with complex amplitude (p ■ 1,2,...,P) 
irifo at mo^ KLP correlator outputs con^esponding to at most L nuittipath signals for each of the /C users 
at each of the P antenna array elements; 

ii) means for weighting the output of each correlator according to the complex conjugate of an estinDate of 
ss the coefficient for a respedive channel and path (h ]^ , ^ c ]^ ;); 

Hi) means for sunriming together the weighted resulting LP complex values for each user, arxl 
Iv) means for forming a /C-vector of the K resulting summations; with 
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b) means for applying a linear combiner to the K-vector to suppress nrultiple acx:es5 interfererKe; and 

c) means for processing the resulting signals to produce estimates of the respective data symbols. 

26. A receiver according to claim 25 wherein said means for the linear combiner con^rises means for forming the prod- 
5 uctof: 

the diagonal anrplrtude matrix (A); 

the equivalent array/channel/code correlation matrix (M); and 
the matrix inverse of the matrix 

10 
IS 

where is the hermitian conjugate of the array/channel/code correlation matrix, 

is the estimated array/channel/code correlation matrix and (a^) is the background noise variance. 

2£ 27. A receiver according to claim 25 wherein the real parts of the K resulting summations are extracted prior to the lin- 
ear combiner. 

28. A receiver comprising: 

30 a bank of correlators for despreading K code-spread data signals S1.S2 Sk bearing respective data 

^i.£»2 each with at most L resolvable muttipath components to form a K-vector, 

means tor applying a linear combiner to the K-vector to suppress muttple access interference; and 

means for processing the output of the linear combiner to produce estimates of the respective data symbols. 

35 29. A receiver according to claim 28 wherein said means for the linear cbmtxner comprises forming the product of: 

the diagonal anplrtude matrix (A); 

the hermitian transpose of the array/channel matrix: arxi the matrix inverse of the matrix 



where 

45 

ft 



so is the complex block-Toepl'rtz correlation matrix o^ is the background noise variance, and I^lp is the identity 

matrix whose dimension is KLP x KLP, 

30. A method of demodulating K direct sequence code-spread data signals subjected to /.-path frequency-selective 
fading, received at an array of P antennas comprising: 

55 

despreading with a bank of fitters matched to L muttipath delayed replicas of each of the K signals received at 
each of the P antennas: 

weigliting the despread signals with the complex conjugate of an estimate of the ray/channel coeff id&it: 
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summing together the weighted, despread signal to form a K-vector; 

applying a linear transform to the K-vector to suppress mutttple access interference; and 

processing the resulting signals to produce estimates of the respective data symbols. 

31 . A method according to claim 30. wherein said linear transform is the product of: 

the cfiagonal amplitude matrix (A); 

the equivalent array/channel/code correlation matrix (M); and 
the matrix inverse of the matrix 

(MA^M'^+iT^M), 

where /W^ is the hermitian conjugate of the array/channel/code correlation matrix. 

(M) 

is the estimated array/channel/code correlation matrix and (o^) is the background noise variance. 

32. A method of demodulating K direct sequence code-spread data signals subjected to L-path frequency-selective 
fading, received ctt an ray of P antennas comprising: 

a) despreading with a bank of filters matched to L multipath delayed replicas of each of the K signals received 
at each of the P antennas to produce a K-vector. 

b) applying a linear transform to the K-vector to combine the nrtatched filter outputs and suppress multiple 
access interfererrce; and 

c) processing the resultant signals to produce estimates of the respective data symbols. 

33. A method according to claim 32 wherein said means tor the linear combiner conprises forming the product of: 

d) the diagonal amplitude matrix (A); 

e) the hermitian transpose of the array/channel matrix (G^; and 

f) the matrix inverse of the matrix 

where 

is the complex block-Toeplitz correlation matrix, is the t>ackground noise variance, and l/^p >s the identity 
matrix whose dimension is KLP x KLP. 

34. A receiver for use with direct-sequence spread spectrum signals, comprising: 

a rake receiver with array processing; arvl 

a muftiuso" detector receiving an output of said rake receiver. 

35. A recaver according to claim 34 wherein said rake receiver with array processing includes: 

a bank of f Qters matched to the nruittipath timing delays of the nuittiple spread spectrum signals; and 
a droflt for weighting the matched filter outputs with the corresponding array/channel estimates. 
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36. A system for detecting the data signals modulated respectively by multiple direct-sequence spread spectrum sig- 
nals, comprising: 

a plurality of space diversity antennas; 
5 a plurality o1 matched filters associated with said arrtennas for providing a plurality of outputs (z^ ... Zp) for at 

least a subset of the multiple signals, a plurality of the multipath components of these signals, and a plurality 
of receiver antennas; 

a circuit for obtaining coherent channel estimates ( h i^ ipC f^jl from the matched filter outputs; 
a second circuit for weighting and combining the matched fitter outputs using the coherent channel estimates; 
10 and 

a third circuit for estimating the respective data symbols for at least a subset of the spread spectrum signals 
white suppressing multiple access interference. 

37. A receiver comprising a processor for executing software, said software being in a computer readable medium and 
16 being for detecting the data signals modulated respectively by multiple direct-sequence spread spectrum signals. 

said software causing said processor to perform at least the steps of: 

a) generating matched fitter outputs for at least a subset of the multiple signals, a plurality of the multf^ath com- 
ponents of these signals, arxf a plurality of receiver antennas; 
20 b) developing coherent channel estimates from the matched f 3ter outputs: 

c] weighting and combining the matched filter outputs using the coherent channel estimates; and 

d) estimating the respective data symbols for at least a subset of the spread spectrum signals while suppress- 
ing multiple access interference. 
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